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ABSTRACT 

From observations collected with the ESPaDOnS spectropolarimeter at the 
Canada- France-Hawaii Telescope (CFHT), we report the detection of Zeeman sig- 
natures on the low-mass classical T Tauri star (cTTS) V2247 Oph. Profile distortions 
and circular polarisation signatures detected in photospheric lines can be interpreted 
as caused by cool spots and magnetic regions at the surface of the star. The large-scale 
field is of moderate strength and highly complex; moreover, both the spot distribution 
and the magnetic field show significant variability on a timescale of only one week, as 
a likely result of strong differential rotation. Both properties make V2247 Oph very 
different from the (more massive) prototypical cTTS BP Tau; we speculate that this 
difference reflects the lower mass of V2247 Oph. 

During our observations, V2247 Oph was in a low-accretion state, with emission 
lines showing only weak levels of circular polarisation; we nevertheless find that ex- 
cess emission apparently concentrates in a mid-latitude region of strong radial field, 
suggesting that it is the footpoint of an accretion funnel. 

The weaker and more complex field that we report on V2247 Oph may share 
similarities with those of very-low-mass late-M dwarfs and potentially explain why low- 
mass cTTSs rotate on average faster than intermediate mass ones. These surprising 
results need confirmation from new independent data sets on V2247 Oph and other 
similar low-mass cTTSs. 
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- stars: individual: V2247 Oph - techniques: spectropolarimetry 
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Whereas our understanding of most phases of stellar evolu- 
tion made considerable progress throughout the twentieth 
century, stellar formation remained rather enigmatic and 
poorly constrained by observations until about two decades 
ago with the advent of the most powerful ground-based and 
space telescopes. One of the major discoveries obtained then 
is that protostellar accretion discs are often associated with 
extremely powerful and highly collimate d jets escap i ng th e 
systems along their rotation axis (e.g., ISnell et al.|[l980l ). 
This finding has revolutionized the field of stellar formation; 
in particular, it provided us with strong evidence that mag- 
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netic fields are playing an important role throughout stellar 
formation. 

Magnetic fields are expected to modify significantly the 
contraction of molecular clou ds into protostellar cores and 
discs (e.g., lAndre et al. I l2008l . for a review). They are also 
expected to influence the next formation stage, at least for 
stars with masses lower than about 2.5 Mq; these protostars 
(called classical T Tauri stars or cTTSs) apparently host 
magnetic fields strong enough to disrupt the central regions 
of their accretion discs, connect the protostars to the inner 
disc rim through discrete accretion funnels and possibly even 
slow down the rotation of protos tars through the res ulting 
star/disc magnetic torque (e.g.. iBouvier et al.ll2007l . for a 
review) . 

While magnetic fields have been repeatedly reported 
over the last decade at the surface of a number of pro- 
totyp ical TTSs (e.g.. IJohns-KruU et al.| [l999: Jo hns-Krulll 
l2007i ). it is only recently that their large-scale topologies 
can be investigated into some det ails through phase-resolved 
spectropolarimetric observations (jPonati et al.ll2007l . l2008bl : 
iHussain et"alll2009l '). This new option offers a direct oppor- 
tunity to pin down the role of magnetic fields in this crucial 
phase of the formation process. In particular, it straightfor- 
wardly allows studying how magnetic topologies of cTTSs 
depend on mass and rotation rate, giving hints on the po- 
tential origin of such fields; moreover, it provides material 
for quantitatively studying how protostars magnetically con- 
nect and interact with their accretion discs. 

Up to now, magnetic topologies of 4 cTTSs have been 
investigated with this method. Among these 4, only BP Tau 
(M* ~ 0.7 Mq) is fully convective and hosts a rather simple, 
dominantly dipolar, large-scale magnetic field, the 3 oth- 
ers being partly convective (as a result of their relatively 
higher mass, in excess of 1.3 Mq) and hosting a more com- 
plex magnetic topology. In particular, this transition from 
simple to complex fields for stars on either sides of the 
full convection limit is reminiscent of t he magnetic proper- 
ties of low-mass ma in-sequence dwarfs (|Morin et al.ll2008bl : 
iDonati et al ]|2008d ): in addition to suggesting that magnetic 
fields of cTTSs are likely of dynamo origin, it also gives po- 
tential hints on why magnetospheric accretion and rotation 
properties of cTTSs are discrepant on both sides of the full 
convection limit (e.g., with higher mass cTTSs rotating more 
quickly in average). 

Extending this magnetic survey to a larger sample of 
cTTSs is needed to go further; this is the exact purpose of 
the MaPP (Magnetic Protostars and Planets) program, in 
the framework of the international MaglcS research initia- 
tive (aimed at qualifying the magnetic properties of stars 
throughout the HR diagram). This survey has been allo- 
cated 690 hr of observing time on the 3.6 m Canada-France- 
Hawaii Telescope (CFHT) over 9 successive semesters (from 
2008b to 2012b). The present paper concentrates on the low- 
masfEl cTTS V2247 Oph (spectral type Ml) in the p Oph 
star forming region (Lynds 1688 dark cloud); with a mass 
of only about 0.35 Mq (see Sec. [2)|, it is less massive than 



^ Throughout this paper, cTTSs with masses lower than 0.5 Mq, 
ranging from 0.5 to 1 Mq and larger than 1 Mq are respectively 
called low-mass, intermediate-mass and high-mass cTTSs. 



all cTTSs yet magnetically imaged and thus appears as an 
ideal candidate to expand our sample towards cooler cTTSs. 



2 V2247 OPH = SR 12 = ROX 21 = HBC 263 

First classified as a TTS with weak emission lines (weak-line 
TTS or wTTS) from its mo derate Ha emission (equivalen t 
width of order 0.4 nm, e.g., iBouvier fc Appenzelieilll992l ). 
V2247 Oph has recently been reported to undergo transient 
episodes of strong er H q emission (equivalent width of up 
to 1.76 nm iLittlefair et a l. 20041 that makes it qualify as 
a cTTS based on that criterion (i.e., with an Hq emission 
larger than 0.9 nm at spectral type Ml, e.g., iFang et al.l 
[2009). It therefore suggests that V2247 Oph is still sur- 
rounded by a (likely low-mass) accretion disc. 

The exact nature and evolutionary stage of this disc is 
however unclear. Given that accretion is still on-going (al- 
though sporadically) , we can assume that the inner disc con- 
tains at least gas. However, V2247 O ph shows no significant 
infrar ed excess up to ab out 10 /im (iGras- Velazquez fc Ravi 
i2005i : [cieiaeiaO |2003) indicating that the disc contains 
very little hot-dust. ISOPHO T measurements indicate large 
excesses at 60 and 100 ^m (|Gras- Velazquez fc Ravi I2005D 
suggesting the presence of an outer disc; however, given 
the large size of the aperture used (120") and the IR- 
bright nebulosity lo cated very close to V2247 Oph (e.g., 
iPadgett et allbOOSi ). this excess may potentially be unre- 
lated to V2247 Oph itself and rather result (at least partly) 
from background/foreground contamination. It thus awaits 
for a confirmation from higher spatial resolution data (e.g., 
Spitzer MIPS observations). 

V2247 Oph is also known to be a binary star w ith a 
companion located at about 0.3" (| Simon et al.lll987l ) from 
the primary star, i.e., at 42 AU assuming a distance of 
139 ± 6 pc for V2247 Oph (e.g.. lMamaie3l2008D . The pres- 
ence of the companion likely implies that the disc features 
a gap at a few tens of AUs, with the inner disc being fed 
by a circumbinary structure (possibly that producing the 
far-IR excesses); this configuration may potentially generate 
lower-than-average and unsteady accretion rates like those 
observed on V2247 Oph. At most of the optical wavelengths 
we are working at in the present paper, the companion is 
outshone by the primary star by more than an order of mag- 
nitude and can therefore be neglected for our analysis. 

Fits to the spec tral energy distribution 
iGras- Velazquez fc Ravll2005h indicate that V2247 Oph has 
a photospheric temperature Teff of 3430 K, a radius _R* 
and a luminosity of 2.2 Rq and 0.61 Lq respectively 
(once corrected for a distance of 139 pc) and suffers a visual 
extinction of Ay — 0.45. This is in rough agreement with 
previous determinations fro m both photometric colors an d 
spectroscopic indices (e.g., IBouvier fc Appenzelie3 Il992l ). 
We thus assume here that T^s ~ 3500 K (with a typical 
error of 150 K) and logL*/LQ — —0.2 (with a t ypical error 
of 0.1 dex). Fitting the evolutionary models of ISiess et al.l 
l|200(]| ) to these parameters, we infer that V2247 Oph is a 
fully convective star with a mass of 0.36 Mq (possible values 
ranging from 0.29 to 0.45 Mq), a radius of 2.0 ± 0.3 Rq 
and an age of about 1 Myr. Note that despite the large 
uncertainty on mass (about 25% and esse ntially refle cting 
the uncertainty on eflective temperature, ISiesd [200lh . we 
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can safely conclude that V2247 Oph is truly a low-mass 
cTTS and in particular that it is significantly less massive 
than the intermediate-mass cTTS BP Tau. In the following, 
we assume that M* = 0.35 M© and Rt, — 2 Kq. 

From photometric modulation presumably caused by 
the presence of cool surface spots going in and out of view 
of Earth-based observers, V2247 Oph is reported to be ro- 
tating with a period of about 3 .5 d ijBouvier fc Bertou3 
ll989l : 1Schevchenko fc Herbstlll998h . Long term photometric 
monitoring moreover reveals that this period is significantly 
chan ging with time (from about 3.4 to 3.6 d. lGrankin et al] 
l2008t ): it su ggests that V2247 Oph (as many other cool active 
stars, e.g.. iHaUl Il99lh is experiencing photospheric shear- 
ing due to latitudinal difi^erential rotation, which generates 
different modulation periods for spots located at different 
latitudes. 

With Rc - Ic = 1.3 (jBouvier fc Appenzellej Il992h . 
V2247 Oph is located at a critical posit ion in an HR di- 
agram of young stars l|Lamm et ahlfSoOsh . For stars bluer 
than this value (i.e. more massive than V2247 Oph), the dis- 
tribution of rotation periods is bimodal, with a main peak 
centred at about 7.5 d and an overall average of about 6.7 d 
in the Orion Nebula Cluster (ONC); for stars redder than 
this threshold, the period distribution is unimodal with an 
average of about 3.3 d in the ONC. With a rotation period 
of 3.5 d, V2247 Oph belongs to the second group while the 
prototypical cTTS BP Tau (whose Rc — Ic and rotation pe- 
riod are respectively equal to 0.9 and 7.6 d) belongs to the 
more slowly-rotating populatior0. 

From the photometric period and projected equato- 
rial rotation velocity (equal to 20.5 ± 0.5 kms~^, see be- 
low, in good agreem ent with the lower accuracy estimate of 
iBouvier et al.l[l986l l. we infer that /?* sini ~ 1.4 Rq, where 
i denotes the inclination of the rotation axis to the line-of- 
sight. We obtain that i = 45±10° ; our imaging study further 
confirms this intermediate value of the inclination angle (see 
SecH. 

Optical veiling, presumably caused by accretion and po- 
tentially weakening the strength of the photospheric spec- 
trum, is rather common for cTTSs and often generates large- 
scale and apparently random variations in the equivalent 
width of photospheric lines. For V2247 Oph however, we see 
only small changes (of order a few percent) in the equiva- 
lent widths of LSD Stokes / profiles, suggesting that veiling 
is minimal; we take this as direct evidence that accretion was 
weak on V2247 Oph during our observations. The spectrum 
of V2247 Oph also exhibits all the usual accretion proxies, 
namely Ha, H/3, the He l D3 line and the emission cores 
of the Ca 11 infrared triplet (IRT) lines. From their equiva- 
lent widths and the corresponding line fluxes (see Sec. [5]), 
and using empirical correlations from the publish ed liter- 
ature (e.g.. iMohantv et al.ll2005l : iFang et al.ll2009l ). we can 
derive a rough estimate of the logarithmic mass accretion 
rate M of V2247 Oph (in Mq yr"^) at the time of our run 
of about -9.8; using the same correlations, we find that M 



^ The exact stellar mass at which this limit occurs somewhat 
depends on the accuracy of both atmospheric and evolutionary 
models; while stud ies suggest it is located at about 0.25 Mq 
llLamm et al.ll2005h , our mass estimate of V2247 Oph would place 
it at about 0.35 M©. 



Table 1. Journal of observations. Columns 1-4 respectively list 
the UT date, the heliocentric Julian date and UT time (both at 
mid-exposure), and the peak signal to noise ratio (per 2.6 kms~^ 
velocity bin) of each observation (i.e., each sequence of 4 X 879 s 
subexposures). Column 5 lists the rms noise level (relative to 
the unpolarized continuum level Ic and per 1.8 kms~^ velocity 
bin) in the circular polarization profile produced by Least-Squares 
Deconvolution (LSD), while column 6 indicates the rotational cy- 
cle associated with each exposure (using the ephemeris given by 
Eq.O. 



Date HJD UT S/N ctlsd Cycle 

(2008) (2,454,000+) (h:m:s) (lO'^'/c) 



Jul 


19 


666.78070 


06:40:26 


110 


9.1 


0.937 


Jul 


20 


667.76238 


06:14:11 


120 


8.5 


1.218 


Jul 


21 


668.77056 


06:26:04 


130 


7.5 


1.506 


Jul 


22 


669.76031 


06:11:25 


140 


6.9 


1.789 


Jul 


22 


669.80591 


07:17:06 


140 


6.9 


1.802 


Jul 


24 


671.75990 


06:11:05 


130 


7.7 


2.360 


Jul 


25 


672.77281 


06:29:47 


130 


7.5 


2.649 


Jul 


26 


673.75570 


06:05:16 


150 


6.3 


2.930 


Jul 


27 


674.75594 


06:05:44 


140 


6.8 


3.216 


Jul 


27 


674.88360 


09:09:35 


130 


7.8 


3.252 


Jul 


28 


675.75706 


06:07:28 


60 


18.7 


3.502 


Jul 


28 


675.84904 


08:19:56 


50 


21.6 


3.528 


Jul 


29 


676.75564 


06:05:32 


120 


8.0 


3.787 


Jul 


30 


677.75556 


06:05:33 


140 


7.0 


4.073 


Jul 


30 


677.87649 


08:59:43 


120 


8.5 


4.108 



can reach 10"'"* Mq yr^^ during s tronger accretion epis odes 
(such as that observed by, e.g., iLittlefair et al.l |2004| ) . As 
expected, this is slightly below the typical values expected 
for low-mass cTTSs (e.g.. iMohantv et alllioosl : iFang et al.l 
120091 ). 



3 OBSERVATIONS 

Spectropolarimetric observations of V2247 Oph were col- 
lected in 2008 July using ESPaDOnS on the CFHT in 
Hawaii. ESPaDOnS collects stellar spectra spanning the 
whole optical domain (from 370 to 1,000 nm) at a resolv- 
ing power of 65,00 (i.e., 4.6 km s~^), in either circular or 
linear polarisation ('Donati"20(II^) . A total of 15 circular po- 
larisation spectra were collected over a period of 12 d; each 
polarisation spectrum consists of 4 individual subexposures 
lasting each 879 s and taken in different polarimeter con- 
figurations to remove all spurious polarisation signatures at 
first order. 

Raw frames are processed with Libre ESpRIT, a fully 
automatic reduction package/pipeline available at CFHT. 
It automatically performs optimal extraction of ESPaDOnS 
unpolarised (Stokes /) and circularly polarise d (Stokes V) 
spectra following the procedure described in iDonati et al.l 
(|l997l ). The velocity step corresponding to CCD pixels is 
about 2.6 kms"'^; however, thanks to the fact that the spec- 
trograph slit is tilted with respect to the CCD lines, spectra 
corresponding to different CCD columns across each order 
feature a different pixel sampling. Libre ESpRIT uses this 
opportunity to carry out optimal extraction of each spec- 
trum on a sampling grid denser than the original CCD sam- 
pling, with a spectral velocity step set to about 0.7 CCD 
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V2247 Oph, LSD profiles, 2008 July 24 




-100 -50 



Velocity (km/s) 

Figure 1. LSD circularly-polarized (Stokes V) and unpolarized 
(Stokes /) profiles of V2247 Oph (top/red, bottom/blue curves 
respectively) on 2008 July 24 (cycle 2.360). The mean polarization 
profile is expanded by a factor of 10 and shifted upwards by 1.08 
for display purposes. 



average noise levels of the resulting LSD signatures are of 
order 8xl0~* per 1.8 kms~^ velocity bin, except for the 2 
spectra taken on Jul 28 whose noise is about twice to thrice 
larger. 

In most spectra, Stokes I LSD profiles show clear dis- 
tortions tracing the presence of cool spots at the surface 
of V2247 Oph; in addition. Stokes V LSD profiles show 
Zeeman signatures with peak-to-peak average amplitudes of 
about 0.5% of the unpolarised continuum (i.e., typically 6 
times larger than the mean noise level). An example of LSD 
profiles is shown in Fig. [T] We note that Stokes I and V 
spectra recorded almost exactly 7 d apart from each other 
(thus corresponding to very similar rotation phases but dif- 
ferent cycles, e.g., cycles 1.218 and 3.216, or cycles 1.789 and 
3.787) show significant differences, larger than those usually 
attributable to small phase shifts (e.g., between cycles 3.216 
and 3.252, or cycles 1.789 and 1.802). We speculate that this 
is likely evidence for strong enough differential rotation at 
the surface of V2247 Oph to change the brightness and mag- 
netic spot patterns at the surface of the star in as little as 
7 d; this option is investigated in Sec. |4] 



pixel (i.e. 1.8 kms~^). All spectra are automatically cor- 
rected of spectral shifts resulting from instrumental effects 
(e.g., mechanical flexures, temperature or pressure varia- 
tions) using telluric lines as a reference. Though not per- 
fect, this procedure provides spectra with a relative radial 
velocity (RV) preci sion of better than 0.030 kms^^ (e.g., 
iDonati et alli2008d ). 

The peak signal-to-noise ratios (per 2.6 kms~^ velocity 
bin) achieved on the collected spectra (i.e., the sequence of 
4 subexposures) range between 110 and 150, with 2 lower- 
quality spectra taken in poor weather conditions. Rotational 
cycles E are computed from heliocentric Julian dates accord- 
ing to the ephemeris: 

HJD = 2454663.5 + 3.5E. (1) 

The full journal of observations is prese nted in Table [l] 

Least-Squares Deconvolution (LSDi lOonati et anil997l ) 
was applied to all observations. The line list we employed 
for LS D is compute d from an Atlas9 LTE model atmo- 
sphere l|KurucJl993l ) and corresponds to a M2 spectral type 
(Teff = 3, 500 K and log g = 3.5) appropriate for V2247 Oph. 
We selected only moderate to strong atomic spectral lines 
whose synthetic profiles had line-to-continuum core depres- 
sions larger than 40% neglecting all non-thermal broadening 
mechanisms. We omitted the spectral regions with strong 
lines formed mostly outside the photosphere, such as the 
Balmer and He lines, and the Ca ll H, K and infrared triplet 
(IRT) lines; we also discarded all spectral windows heavily 
crowded with telluric lines. Molecular lines are not included 
in the list, their wavelengths, intensities and Zeeman sensi- 
tivities being not empirically computable with enough ac- 
curacy; we assume that their indirect contribution to LSD 
profiles is not concentrated at discrete velocities but sta- 
tistically distributed over a wide domain, thus inducing no 
significant change nor variability in the shape of Stokes / 
and V LSD profiles. Altogether, about 4,800 spectral fea- 
tures are used in this process, most of them from Fe I. Ex- 
pressed in units of the unpolarised continuum level 7c, the 



4 BRIGHTNESS AND MAGNETIC IMAGING 
4.1 Imaging method 

To model the cool spots and magnetic topology at the 
surf ace of V 2 247 O ph, we use the new imaging code of 
Do nati et al.l (|2006bh where the field is described through 
spherical-harmonics expansions and the brightness distribu- 
tion as a series of independent pixels. In its newest imple- 
mentation, our code can also simultaneously recover the dis- 
tribution of accretion spots as derived from the profiles of 
emission lines formed at the base of accretion funnels (e.g., 
the narrow emission components of the IRT lines). Since 
V2247 Oph was experiencing only weak accretion at the 
time of our observations, we keep this part of the modelling 
separate and describe it in the next section. 

Our imaging code uses the principles of maximum en- 
tropy to retrieve the simplest magnetic image and brightness 
distribution compatible with the series of rotationally modu- 
lated Stokes I and V LSD profiles. More specifically, the field 
is divided into its radial-poloidal, non-radial-poloidal and 
toroidal components, each of them described as a spherical- 
harmonics expansion (e.g.. iDonati et al.|[2006b[ ): given the 
different rotational modulation of Zeeman signatures that 
poloidal and toroidal fields generate, our imaging code ap- 
pears particularly useful and efficient at producing dynamo- 
relevant diagnostics about the large-scale magnetic topolo- 
gies at the surface of late-type stars. The brightness distri- 
bution is simply modeled as a set of independent pixels de- 
scribing the contrast with respect to the quiet photosphere; 
for this particular analysis, we only allow the code to recover 
cool spots (i.e., darker than the photosphere). 

The reconstruction process is iterative and proceeds by 
comparing at each step the synthetic profiles corresponding 
to the current image with those of the observed data set. 
To compute the synthetic Stokes I and V profiles, we divide 
the surface of the star into small grid cells (typically a few 
thousands), work out the specific contribution of each grid 
cell to the Stokes profiles (given the brightness, magnetic 
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field strength and orientation within each grid cell, as well 
as the cell RV, location and projected area) and finally sum 
up contributions of all cells. Since the problem is partly ill- 
posed, we stabilise the inversion process by using an entropy 
criterion (applied to the spherical harmonics coefficients and 
to the brightness image pixels) aimed at selecting the image 
with minimum information among all those compatible with 
the data. 

To describe the local profiles, we 

use Unno- Rachkovsky ' s eq uations (e.g., 

iLandi degl'Innocenti &: Landolfil |2004| ) known to pro- 
vide a fairly accurate description of Stokes / and V profiles 
in the presence of magnetic fields. We set the central wave- 
length, Doppler width and Lande factor of our equivalent 
line to 700 nm, 3.4 kms^^ and 1.2 respectively and adjust 
the average line equivalent width to the observed value. 
We then derive both the rotational broadening vsini and 
the radial velocity Urad of V2247 Oph by fitting our series 
of LSD Stokes I profiles and selecting the values that 
minimise the image information content (for a given quality 
of the fit to the data). We find vsini = 20.5 ± 0.5 kms~^ 
and Vrad = —6.3 ± 0.5 kms~^. We can also obtain a rough 
estimate of the inclination angle by selecting the one 
producing images with minimum information content; the 
value we derive (when fitting both Stokes I and V data 
simultaneously) is 45 ± 10°, in good agreement with our 
initial estimate (see Sec.[2J. 

As noted above, the LSD profiles (and in particular 
their apparent deviation from strict rotational modulation) 
suggest that the brightness and magnetic distributions may 
have changed significantly over the course of our run. By 
comparison with other active cool stars, it is unlikely that 
the field would vary significantly over a period of one week; 
the probable cause for this change is thus surface differen- 
tial rotation. This conclusion is independently confirmed by 
the significant changes in photometri c period reported fo r 
V2247 Oph over the last 2 decades l|Grankin et al.ll2008h . 
To incorpora te differential rotation i nto the mode l ling, w e 
proceed as in lDonati et al.1 l|2003t ) and lMorin et al.1 (l2008al '). 
i.e., by assuming that the rotation rate at the surface of the 
star is varying with latitude 9 as f2cq — dQ sin^ 6 where Qcq is 
the rotation rate at the equator dfl the difference in rotation 
rate between the equator and the pole. When computing the 
synthetic profiles, we use this law to work out the longitude 
shift of each cell at each observing epoch with respect to its 
location at the median observing epoch at which the field 
is reconstructed (at rotation cycle 2 in the present case and 
in the ephemeris of Eq.[Tl i.e., at HJD 2454670.5); this way, 
we can correctly evaluate the true spectral contributions of 
all cells at all epochs. 

For each pair of Qcq and dQ values within a range of ac- 
ceptable values, we then derive, from the complete data set, 
the corresponding magnetic topology (at a given information 
content) and the associated reduced chi-squared Xi a-t which 
modelled spectra fit observations. By fi tting a paraboloid t o 
the Xr surface derived in this process (iDonati et al.ll200^ . 
we can easily infer the magnetic topology that yields the 
best fit to the data along with the corresponding differential 
rotation parameters and error bars. This process has proved 
reliable for esti mating surface diffe rential rotation on mag- 
netic stars (e.g.. iDonati et aLll2003l ). 



4.2 Results 

The first result is that both Stokes I and V LSD pro- 
files confirm that differential rotation is indeed significantly 
shearing the photosphere of V2247 Oph. In both cases, the 
maps describing how Xr varies as a function of f2cq and dQ 
(at constant information content) show a clear paraboloid 
and a well defined minimum (see Fig. [2]). The differential 
rotation parameters producing an optimal fit to the data 
are respectively equal to fioq = 1.89 ± 0.02 radd~^ and 
dQ = 0.34 ± 0.04 radd"^ for Stokes / data and Qcq = 
1.85±0.03 radd"^ and dQ = 0.32±0.05 radd"^ for Stokes V 
data. Both results are compatible within about one a of their 
weighted average {Qaq ~ 1-87 radd'^, dfl — 0.33 radd~^, 
see middle panel of Fig. [2)) that we take as reference here- 
after. The corresponding rotation periods for the equator 
and pole are 3.36 and 4.08 d respectively, bracketing the 
repo rted photo metric periods ranging from 3.4 to 3.6 d 
( Gr ankin et al. 2008); in this context, the period with which 
we phased our data (3.5 d) corresponds to a latitude of about 
30°. 

We carried out several tests to confirm that differen- 
tial rotation is indeed detected. For instance, trying to fit 
Stokes V LSD profiles to a unit reduced chi-squared Xi as- 
suming solid body rotation is hardly possible. The smallest 
Xi we get is about 13% higher (i.e., x^ larger by 70 for 
540 data points) than the minimum obtained when assum- 
ing differential rotation at equal information content, and 
the corresponding rotation period (3.75 d) is significantly 
larger than the average photometric period; Xv even gets up 
to 25% higher (i.e., x^ larger by 135) than the minimum 
when imposing solid body rotation at a period of 3.5 d. We 
also investigated whether the observed differential rotation 
could actually be a spurious signal resulting from intrinsic 
variability (other than differential rotation); note that this 
option is unlikely as there is little chance for intrinsic vari- 
ability to mimic a fast equator and a slow pole in the same 
way in 2 different and fully-independent data sets. We nev- 
ertheless checked it by splitting our data into 2 subsets and 
reconstructing individual images for each subset (assuming 
solid body rotation); although this method has limited di- 
agnostic power (the fewer the spectra in one subset the less 
accurate the recovered image), it is nevertheless clear that 
both images include similar features at slightly different po- 
sitions and different contrasts. We therefore believe that our 
detection is real. 

The surface brightness and magnetic map we derive 
(see Fig. ^ corresponds to a unit Xr fit to the data (see 
Fig. |4)| , the initial Xi (corresponding to an unspotted and 
non-magnetic star) being about 3.4 for 1080 data points in 
total. The brightness map is fairly simple and includes 2 
main features at phases 0.30 and 0.60 with latitudes rang- 
ing from 30 to 70°. In particular, V2247 Oph does not fea- 
ture a polar cap like the conspicu ous ones found, e.g ., on 
the high er mass cTTS V2129 Oph (|Donati et al.ll2007l ) and 
BP Tau l|Donati et al.ll2008bl ). Spots large and dark enough 
to affect significantly the shape of Stokes / LSD profiles do 
not cover more than 3.5% of the stellar surface, at least at 
the time of our observations. Given the significant rotational 
broadening in the spectrum of V2247 Oph, this result is well 
constrained and reliable and cannot be attributed to poten- 
tial defects of surface imaging. 
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Figure 2. Variation of x? as a function of f2cq and dO,, derived from the modelling of our Stokes / (left) and Stokes V (right) data at 
constant information content; the middle plot shows the corresponding confidence ellipses (with Icr and 3cr ellipses shown in full and 
dashed lines). For the left and right plots, the outer contour corresponds to a 2.2% increase in Xv corresponding to a 3(7 ellipse for both 
parameters as a pair. 
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Figure 3. Surface brightness (top left) and magnetic topology of V2247 Oph in 2008 July (HJD 2454670.5), simultaneously derived 
from the Stokes / & V data shown in Fig. U] The radial, azimuthal and meridional components of the field are displayed from top right 
to bottom left (with magnetic flux values labelled in G). The star is shown in flattened polar projection down to latitudes of —30°, with 
the equator depicted as a bold circle and parallels as dashed circles. Radial ticks around each plot indicate phases of observations. 
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Figure 4. Maximum-entropy fit (thin red line) to the observed Stokes / (left) and Stoices V (right) LSD photospheric profiles (thick 
black line) of V2247 Oph. The rotation cycle of each observation (as listed in Table 1) and la error bars (for Stokes V profiles) are also 
shown next to each profile. 



The large-scale magnetic field of V2247 Oph that we 
recover i s also fairly d ifferent from those of V2129 Oph and 
BP Tau tPonati et al.| [2007. 2008tii); not only is the detected 
magnetic flux much weaker (about 300 G in average) but 
its topology is also fairly different, with a strong toroidal 
component (totalling about 60% of the magnetic energy, as 
obvious from Fig.[3| and a predominantly non-axisymmetric 
poloidal component (less than 30% of the poloidal field en- 
ergy concentrating in modes with m < £/2). Moreover, the 
poloidal field is complex with about half of the poloidal field 
energy concentrating in modes with £ > 2; this is fairly ob- 
vious from the radial field component, featuring nearby re- 
gions of opposite polarities (see Fig. [Sjl . 

To investigate further the large-scale magnetic topology 
of V2247 Oph, we forced the code to recover fields that are 
either mainly symmetric or antisymmetric with respect to 
the centre of the star; this is achieved by favouring spherical 
harmonics modes of even or odd I values respectively. We 
find that a symmetric field configuration is marginally more 
likely; in both cases though, the data can easily be fitted 
down to noise level and the recovered topologies in the vis- 
ible hemisphere are very similar. The dipole component of 
the poloidal field is weak and strongly tilted with respect to 
the rotation axis; its polar strength ranges from 80 to 100 G 
and concentrates about 25-30% of the poloidal field energy 
(depending on the assumed symmetry properties). 

Note that the intrinsic weakness of the dipole compo- 
nent on V2247 Oph is a reliable and well constrained result; 



in particular, it cannot be attributed to observational prob- 
lems, e.g., to magnetic flux hiding in large high-contrast dark 
spots (absent on V2247 Oph at the time of our observations) 
or to flux cancellation from regions of opposite polarities 
(easily detected in the radial held component). 



5 ACTIVITY & ACCRETION 

As usual for cTTSs, V2247 Oph exhibits in its optical spec- 
trum a number of emission lines known to be proxies of 
activity and/or accretion - and in particular, the Balmer 
lines Hq and H/3, the He I D3 line and the 849.8 nm line 
of the Ca 11 IRT. The average equivalent widths that we 
measure in these proxies are equal to 137, 140, 8.5 and 
12 kms~^ respectively, or equivalently 0.30, 0.23, 0.017 and 
0.035 nm; they correspond to logarithmic line fluxes of 6.2, 
5.8, 4.9 and 5.4 (with fluxes in ergs~^cm~^), or equiva- 
lently -4.6, -5.0, -5.9 and -5.4 (with fluxes in L©). From 
these, we derive an average logarithmic accretion luminos- 
ity (in L0) of —3.2 ± 0.3 and a logarithmic mass accretion 
rate (in Mq yr~^) of -9.8 ± 0.3 jPang et alj|2009l ). 

In this paper, we concentrate mostly on 3 proxies. Ha, 
H/3 and the LSD average of all 3 Ca 11 IRT lines (averaging 
out to an equivalent width of about 14 kms~^ or 0.04 nm). 
The He I D3 line is dim and noisy and hardly usable for the 
present analysis. The complete series of proflles for all three 
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proxies are shown in Fig. (5] small Stokes V signatures are 
also detected in conjunction with the IRT lines. 

As obvious from Fig. [S] the accretion proxies feature 
a significant level of intrinsic variability. For instance, the 
amount of emission at cycle 1.789 is clearly different from 
that 2 full rotation cycles later (at cycle 3.787) in all 3 prox- 
ies; the emission peak detected at cycles 2.360 and 2.649 
also seems fairly difficult to reconcile with emission fiuxes 
at nearby phases in both the previous and following rota- 
tion cycles. To separate intrinsic variability from rotational 
modulation (generating mostly low frequency periodic varia- 
tions) , we propose the following procedure. We start by com- 
puting the equivalent widths of all 3 proxies as a function of 
rotation cycle (excluding the low quality data collected on 
28 July) and fit them with a low-frequency periodic wave - 
hopefully retaining rotational modulation only and filtering 
out intrinsic variability; we then simply rescale all profiles to 
their fitted equivalent widths. While obviously no more than 
approximate, this procedure has the advantage of being very 
straightforward yet reasonably efficienl|f|. 

The periodic wave we use to fit the equivalent widths is 
of the form: 

ai cos{2iTvE/p) + h sm{2iTxE/p) + c (2) 
where E denotes the rotation cycle of each data point 
(as in Eq. [T]) and p, ai, bi and c are the 6 free parameters 
of our model wave. We find that for the 3 proxies, the fit 
residuals show a very clear minimum for p ~ 1.08, i.e., for a 
period of about 3.8 d. Assuming that the differential rota- 
tion law derived in Sec. |4] also applies in the chromosphere 
(where accretion proxies supposedly form), it suggests that 
active and/or accreting regions mainly concentrate at inter- 
mediate latitudes of about 50° . An example of the modelled 
rotational modulation in the particular case of the IRT is 
shown in Fig. |6] 

By rescaling all IRT profiles to their fitted equivalent 
widths, we obtain a data set that can be added to that 
from photospheric lines (see Sec. |4} and further constrain 
the modelling of V2 247 Oph. In that respect , we proceed 
as m previous work (|Donati et al.ll2007l . l2008bl ) and assume 
that IRT emission includes both a quiet chromospheric con- 
tribution from the whole stellar surface as well as a stronger 
but more localised contribution from accretion regions. In 
previous work, we assumed (as a first guess) that accretion 
spots coincide with dark photospheric spots; we no longer 
use this approximation here and recover 2 completely inde- 
pendent distributions (in addition to the magnetic field) , one 
giving the photospheric surface brightness (as in Sec.|4]) and 
another one describing the amount of excess emission from 
accretion spots. The local emission line (modelled again us- 
ing Unno-Rachkovsky's equations) features a central wave- 
length, Doppler width and Lande factor of 850 nm, 7 kms~^ 
and 1.0 respectively. The intrinsic emission within accretion 
spots is (arbitrarily) assumed to be 10 times that within 
the quiet photosphere; the exact value of this ratio has lit- 
tle impact on the result, e.g., with smaller ratios producing 

^ Our scaling procedure has little impact on the reconstructed 
maps and mostly aims at filtering out the discrepancies between 
the synthetic and observed profiles that the model was not de- 
signed to reproduce; in other words, it essentially makes the code 
convergence quicker and more reliable by emphasising the rele- 
vant contrasts in the landscape. 



only more contrasted accretion regions without noticeably 
changing their distribution. 

The surface distribution of excess IRT emission we ob- 
tain is shown in Fig. [7] while the corresponding fit to the 
Stokes I and V profiles is shown in Fig. [S] We find that ex- 
cess IRT emission is concentrated in one main region located 
at phase 0.45 and latitude 45°, and covering about 1.5% of 
the stellar surface; unsurprisingly, this result is compatible 
with the observed modulation of the Ca ll IRT core emission 
(peaking at phase 0.45 and fluctuating with a period corre- 
sponding to intermediate latitudes, see Fig. |6} . The bright- 
ness and magnetic maps at photospheric level (as well as 
the corresponding fits to Stokes / and V LSD profiles) are 
virtually identical to those already shown in Sec. [l] The 
shapes of Stokes I IRT profiles (and in particular the asym- 
metries) are well reproduced in our model, ensuring at the 
same time that the corresponding RVs (observed to vary by 
about ±2 kms~^ about the line centre) are also grossly fit- 
ted (see Fig. |9]). Given how well correlated the equivalent 
width and RV curves are for the IRT line. Ha and H/3, we 
assume that the excess IRT emission spot we detect is also 
what generates the observed excess Ha and H/3 emission. 

Worth noticing is that the average excess-emission spot 
latitude I that we obtain (about 45°) is much less than that 
one would derive from naively assuming that all emission 
comes from the accretion spot and hence that the semi- 
amplitude of the line RV variations is roughly given by 
wsinicosZ (which would imply / ~ 85°). The excess emis- 
sion from accretion spots (generating the RV variations) is 
actually mixed with the (constant and dominant) contribu- 
tion from the quiet chromosphere (the light-blue line in the 
left panel of Fig. [8]), thereby strongly reducing the observed 
amplitude of the RV variations. The shape of IRT profiles 
visually demonstrates this; the IRT line is not wiggling back 
and forth in velocity space as a whole (as expected from 
emission concentrating in a localised accretion spot) but fea- 
tures a roughly constant centroid and width at about 10% 
height (as expected when the whole star significantly con- 
tributes to emission through a quiet chromosphere) . We thus 
(unsurprisingly) conclude that, fitting emission line profiles 
(whenever possible) is much more efficient and reliable than 
modelling RV variations only, for mapping accretion spots 
at the surfaces of cTTSs. 



6 SUMMARY AND DISCUSSION 

We presented in this paper the first magnetic map for the 
low-mass cTTS V2247 Oph. This is the lowest mass cTTS 
to be magnetically imaged yet; despite it being in a low- 
state of accretion (with logM ~ —9.8 with M expressed 
in Mq yr~^), this study brings a number of new and unex- 
pected results. 

The brightness distribution we derive at the surface of 
V2247 Oph (from the rotational modulation of Stokes / 
photospheric LSD profiles) is relatively simple, with dark 
spots covering about 3.5% of the total surface; in particu- 
lar, it does not feature a high ly-contrasted dark polar spot 
as seen, e.g., on V2 129 Oph jDonati et al.ll2007h. BP Ta u 



(|Donati et al.ll2008bD or other cTTSs (|Hussain etaLll2009l ). 

It is however roughly similar to what is found on low-mass 
fuUy-convective dwarfs where large-scale spots are usually 
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Figure 5. Variation of Ha (left), H/3 (middle) and IRT (right) profiles throughout our observing run. In each case, the average profile 
over the run is shown in red. The rotation cycle of each observation (as listed in Table I) are also shown next to each profile. 




Rotation cycle 

Figure 6. Measured (red open circles) and fitted (blue crosses) 
equivalent widths of the Ca II IRT LSD profiles throughout our 
observing run. The model wave (green line) providing the best fit 
to the data features a period of about 3.8 d. 
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Figure 7. Distribution of excess Ca IRT emission at the surface 
of V2247 Oph. 



low-contrast and cover only a s mall fraction of the stellar 
surface (e.g.. iMorin et al.ll2008al '). 

We also detect clear Zeeman signatures from 
V2247 Oph, tracing a rather complex and moderately 
strong multipolar large-scale magnetic topology with an av- 
erage strength of 300 G and a mixed amount of poloidal and 
toroidal field; the poloidal field is mostly non-axisymmetric 
and features a weak (tilted) dipole component of about 
80-100 G. Again, this is fairly different from what is seen 
on intermediate-mass cTTSs like BP Tau, where the field 
is far simpler, predominantly poloid al, axisymmetric wit h 
a dipole component exceeding 1 kG ijPonati et al.l[2008bl ). 
This is also in contrast to low-m ass fuUy-convective mid-M 
dwarfs (e.g.. iMorin et al.ll2008bl ) that almost always exhibit 



very intense, mainly potential large-scale fields with a 
simple topology roughly aligned with the rotation axis. 
New results however indicate that a significant fraction of 
very-low-mass dwarfs (late-M dwarfs with M* < 0.2 M©, 
i.e., below about half the mass threshold for full convection 
on the main sequence) are found to host moderate and com- 
plex non-axisymmetric large-scale magnetic fields (Morin 
et al 2009, in preparation), very different from the strong 
and simple fields of mid-M dwarfs. Although V2247 Oph is 
significantly more massive (at 0.35 Mq) than late-M dwarfs, 
the difference in evolutionary stage may easily compensate 
for the discrepancy; the mass of V2247 Oph is indeed also 
lower than half the mass threshold for full convection at 
an age of about 1 Myr. Further observations are of course 
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Figure 8. Maximum-entropy fit (thin red line) to the observed Stokes / (left) and Stokes V (right) IRT emission profiles (thick black 
line) of V2247 Oph; the light-blue curve shows the contribution of the quiet chromosphere. Note that all profiles were rescaled following 
the procedure described in the text. The rotation cycle of each observation and Icr error bars (for Stokes V profiles) are also shown next 
to each profile. 
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Rotation cycle 

Figure 9. Measured RVs of the IRT (red circles), Ha (green 
crosses) and H/3 (blue squares) profiles throughout our observ- 
ing run. The light-blue line shows the RV curve that our model 
predicts. 



needed to confirm whether this analogy truly holds; if so, 
it would bring additional evidence that magnetic fields of 
intermediate- and low-mass cTTSs are similar in nature to 
those of mid-M and late-M dwarfs and are thus produced 
through dynamo processes (rather than being fossil leftovers 
from an earlier formation stage). 



In addition, we find that the spot distribution and 
the magnetic field of V2247 Oph evolve on a very short 
timescale, of order one week; in both cases, this evolution 
is compatible with surface difi'erential rotation shearing the 
photosphere 5 to 6 times faster than on the Sun, with the 
equator lapping the pole by one complete rotation every 
19 ± 2 d. This is again fairly unexpected since all fuUy- 
convective stars have been reported to show little to no 
differentia l rotation up to now (e.g., iDonati et al.l l2006al : 
iMorin et~a l. 2008a, b). Only late F stars with very shallow 
convective zones are yet known to exhibit a si milar degree of 
photospheric shear (e.g.. lOonati et al]|2008al ). Further con- 
firmation from new data sets on V2247 Oph (and other sim- 
ilar low-mass cTTSs) are thus needed to validate this sur- 
prising result. 

The amount of differential rotation we estimate is com- 
patible with the range of photometric periods measured on 
V2247 Oph over the la st 2 decades (from 3.4 to 3.6 d, imply- 
ing dVL > 0.1 radd~\ [Crankin et aLlbOOSl ). Given the dif- 
ferential rotation law we derive, these periods suggest that 
spots preferentially cluster at low to intermediate latitudes 
(15-40°) on V2247 Oph; this is at least qualitatively compat- 
ible with the fact that no high-contrast polar spot is present 
at the surface. This is the first time that differential rotation 
is detected on a cTTS; the recent claim (|Herbst et al.ll2006l ) 
that the high-mass cTTS HBC 338 (spectral type G9) is dif- 
ferentially rotating (given the observed variations of the pho- 
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tometric period) needs confirmation, cTTSs with massive 
accretion discs being prone to photometric perturbations 
(including sudden changes of the hght-curve period) likely 
caused by accretion (e.g., ISimon et al.lll990l : iDonati et al.l 
l2008bl ) rather than resulting from differential rotation. 

Finally, we find that V2247 Oph features a region of 
excess IRT (and Ha and H/3) emission at phase 0.45 and 
intermediate latitude. There is no obvious correlation be- 
tween the location of this emission region and the cool spots 
mapped from the distorted LSD profiles of photospheric 
lines; however, we find that it roughly coincides with one 
of the strongest field region detected on V2247 Oph, the 
positive radial field spot whose magnetic fiux reaches about 
0.5 kG. We speculate that this spot may be the footpoint of 
an accretion funnel linking the star to its accretion disc at 
the time of our observations. If confirmed, we expect such 
accretion funnels to be rather short lived, with differential 
rotation continuously distorting the magnetic connections 
between the star and the disc as well as the stellar field itself. 
In particular, this effect may be partly responsible for the 
sporadic accretion episodes observed on V2247 Oph, with 
more intense accretion episodes occurring when the mag- 
netic topology linking the star to the disc is favouring ac- 
cretion. 

In addition to the differences in their magnetic topolo- 
gies, V2247 Oph and BP Tau show quite distinct rotational 
properties; as emphasised in Sec. [2] ea ch of them belong s 
to a distinct population of young stars (|Lamm et al.|[2005f ), 
V2247 Oph to the low-mass cTTSs with short rotation pe- 
riod (of about half a week in average) while BP Tau to 
the intermediate-mass cTTSs with long rotation periods (of 
about one week in average) . These two issues may relate to- 
gether, i.e., the magnetic topology may have an impact on 
the rotation period through a coupling mechanism such as, 
e.g., the wel l-known disc-locking picture (Camonzind 1 990l : 
lKonigllll99lh . In this model, the rotation period at the sur- 
face of the star settles close to the Keplerian period of the 
point in the disc (called the Alfven radius and noted ta) 
where the ram pressure of the average accretion fiow roughly 
equals the magnetic pressure of the largest-scale (i.e., th e 
dipole compo nent) poloidal f i eld (e .g., lBessolaz et"al]l2008f ). 
Using Eq. 1 of lBessolaz et al.l (|2008l ) (with _B* set to the aver- 
age dipole field strength at the rotation equator, i.e., about 
40 G for V2247 Oph), we estimate ta ~ 3.9 i?*; we can 
also work out that the radius at which the Keplerian period 
matches the average rotation period (called the corotation 
radius and noted rc) is equal to rc — 3.4 7?* for V2247 Oph. 
Both estimates roughly agree with each other, suggesting a 
reasonable agreement with the disc-locking model predic- 
tions. 

Comparing directly with BP Tau however (and assum- 
ing a logarithmic m ass accretion rate o f -7.5 taken from 
the literature as in iDonati et all l20Q8bt ). the model pre- 
dicts that both stars should have similar va and thus be- 
have similarly with respect to the disc-locking mechanism, 
the 15-fold weaker dipole field of V2247 Oph being almost 
compensated by the 200-fold weaker mass accretion rate; re- 
estimating log A/ for BP Tau using the same procedure as 
for V2247 Oph, we actually find a much smaller value, equal 
to —8.6 ± 0.4. While it is not clear whether this new esti- 
mate is more accurate than the one used previously, it is at 
least more consistent with the one we derived for V2247 Oph 



when it comes to comparing both stars; in this case, we ob- 
tain that rA ~ 7.1 for BP Tau (taking ~ 600 G) in 
reasonable agreement with rc — 7.4 

Obviously, more stars of different masses and rotation 
periods, and with consistent estimates of _B* and M, are 
needed to assess the validity of the disc-locking picture; if 
confirmed, it may ultimately explain why low-mass cTTSs 
rotate in average faster than their intermediate mass equiv- 
alents. The MaPP program should bring new results along 
this line in the near future. 
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